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First-principles  electronic  structure  calculations  have  been  very  useful  in  understanding  some 
of  the  properties  of  the  new  iron-based  superconductors.  Further  explorations  of  the  role  of  the 
individual  atomic  orbitals  in  explaining  various  aspects  of  research  in  these  materials,  including 
experimental  work,  would  benefit  from  the  availability  of  a  tight-binding(TB)  Hamiltonian  that 
reproduces  accurately  the  first-principles  band  structure  results.  In  this  work  we  have  used  the 
NRL-TB  method  to  construct  a  TB  Hamiltonian  from  Linearized  Augmented  Plane  Wave(LAPW) 
results.  Our  TB  model  includes  the  Fe  d-orbitals,  and  the  p-orbitals  from  both  As  and  O  for  the 
prototype  material  LaOFeAs.  The  resulting  TB  band  structure  agrees  well  with  that  of  the  LAPW 
calculations  in  from  2.7  eV  below  to  0.8  eV  above  the  Fermi  level,  ef,  and  the  Fermi  surface  matches 
perfectly  to  that  of  the  LAPW.  The  TB  densities  of  states(DOS)  are  also  in  very  good  agreement 
with  those  from  the  LAPW  in  the  above  energy  range,  including  the  per  orbital  decomposition. 
We  use  our  results  to  provide  insights  on  the  existence  of  a  pseudogap  in  the  DOS  just  above  the 
Fermi  level.  We  have  also  performed  a  separate  TB  fit  to  a  database  of  LAPW  results  as  a  function 
of  volume  and  with  variations  of  the  As  positions.  This  fit  although  less  accurate  regarding  the 
band  structure  near  ef,  reproduces  the  LAPW  total  energies  very  well  and  has  transferability  to 
non-fitted  energies. 


I.  INTRODUCTION 

The  discovery  of  superconductivity  with  critical  tem¬ 
perature  Tc  of  about  55K  in  the  iron  compounds  named 
iron-Oxypnictidesii~-  has  brought  to  the  field  excitement 
comparable  to  that  created  by  the  high  Tc  cuprates. 
The  prototype  formula  for  these  materials  is  LaOFeAs 
with  two  distinct  layers  of  LaO  and  FeAs.  The  room- 
temperature  crystal  structure  is  tetragonal  and  under¬ 
goes  a  structural  distortion  to  orthorhombic  at  low  tem¬ 
peratures.  The  transition  temperature  is  modulated  by 
electron  doping  (F  substituting  for  O)  and  hole  doping 
(Sr  substituting  for  La).  Other  substitutions  may  occur, 
e.g.,  replacing  Fe  by  Ni  or  As  by  P. 

Many  experimental  and  theoretical  investigations  have 
been  performed  without  reaching  a  consensus  as  to 
whether  these  materials  are  similar  to  the  cuprates,  with 
their  superconducting  mechanism  not  yet  understood,  or 
if  the  essential  physics  is  different,  indicating  that  this 
is  yet  another  class  of  superconductors.  Questions  to 
be  answered  include  the  nature  of  the  normal  state,  the 
symmetry  of  the  superconducting  state  and,  of  course, 
the  origin  of  the  pairing  interaction. 

Density  functional  theory  (DFT)  calculations^  have 
been  at  the  center  of  the  theoretical  investigations.  For 
LaOFeAs,  Singh  and  Du£  found  a  high  density  of  states 
at  the  Fermi  level,  N(£^),  and  a  low  carrier  concentra¬ 
tion.  We  note  that  this  is  different  from  the  cuprates 
which  have  a  low  N(£i?)  and  while  displaying  low  car¬ 
rier  concentrations  they  are  characterized  by  a  half-filled 


band  near  ef ■  On  the  other  hand,  Singh  and  Du  found 
a  high  N(£jQ  with  antiferromagnetic  fluctuations,  which 
has  a  definite  similarity  with  the  cuprates.  On  the  same 
theme,  competing  antiferromagnetism  and  superconduc¬ 
tivity  in  the  doped  system  is  suggested  by  Yildirim^  as 
breaking  the  tetragonal  symmetry  causing  a  structural 
distortion. 

Since  standard  DFT  is  very  expensive  computationally 
there  is  a  need  for  developing  tight-binding  (TB)  mod¬ 
els  that  can  be  the  starting  point  for  carrying  out  fur¬ 
ther  investigations  using  many-body  techniques,  such  as 
multiband  Hubbard  models.  A  few  TB  approaches  have 
appeared  in  the  literature,  such  as  the  study  of  Kuroki  et 
al.7',  which  is  based  on  an  Fe-only  d-band  Hamiltonian. 
These  authors  applied  the  random  phase  approximation 
to  obtain  spin  and  charge  susceptibilities,  concluding  that 
an  unconventional  s-wave  pairing  is  in  play. 

Having  examined  the  details  of  the  DFT  calculations 
for  LaOFeAs,  we  identified  the  following  features  of  the 
energy  bands  emerging  from  a  wave  function  analysis. 
Starting  from  the  lower  bands  we  find  O  p-states  that,  for 
higher  energies,  hybridize  with  As  p-states.  Hybridiza¬ 
tion  with  Fe  d-states  occurs  at  the  top  of  the  As  p-states 
and  then  hybridization  with  La  states  appears  high  above 
the  Fermi  level.  Therefore,  it  becomes  clear  that  TB 
models  that  ignore  the  other  elements  and  use  only  the 
Fe  d-orbitals  are  not  reproducing  the  band  structure  of 
LaOFeAs  accurately  enough. 

Furthermore,  a  recent  paper  by  Manousakis  et  al£- 
builds  a  TB  Hamiltonian  fitted  to  DFT  results  using  in 
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addition  to  the  Fe  d-orbitals  p  orbitals  of  As.  Based 
on  their  TB  Hamiltonian  these  authors  report  that  the 
effective  Hamiltonian,  in  the  strong  on-site  Coulomb- 
repulsion  limit,  operates  on  three  distinct  subspaces  cou¬ 
pled  through  Hund’s  rule.  They  also  argue  that  the 
observed  spin-density-wave  order  minimizes  the  ground 
state  energy  of  the  Hamiltonian.  These  conclusions  could 
be  a  correct  speculation  of  the  physics  in  this  mate¬ 
rial.  However,  their  calculations  are  based  on  a  TB 
Hamiltonian  that  is  not  accurately  derived  from  the  first- 
principles  data.  Although  they  performed  their  TB  fit 
only  near  the  Fermi  level,  their  results  do  not  reproduce 
the  energy  bands  well  enough  as  it  can  be  seen  in  their 
Fig  4.  Therefore,  the  conclusions  of  this  paper  based  on  a 
poorly  constructed  TB  Hamiltonian  can  only  be  consid¬ 
ered  as  a  speculation.  More  recently  a  paper  by  Eschrig 
and  Koepernik—  presents  a  minimal  basis  TB  Hamilto¬ 
nian  for  LaOFeAs,  as  well  as  for  other  structure  types 
of  the  Fe  superconductors,  that  is  based  on  an  elegant 
TB  theory.  However,  this  Hamiltonian  has  quantitative 
agreement  to  first-principles  results  only  in  a  smaller  win¬ 
dow  around  E^  than  the  fit  we  present  here. 

In  this  work  we  have  used  the  Naval  Research  Labo¬ 
ratory  Tight-Binding  (NRL-TB)10>ii  method  to  fit  our 
linearized  augmented  plane  wave  (LAPW)  resultsi2~-iito 
a  TB  basis  with  the  aim  of  reproducing  the  band  struc¬ 
ture  very  accurately.  We  have  included  the  d  orbitals  of 
Fe,  the  p  orbitals  of  As,  and  the  p  orbitals  of  O,  leaving 
out  the  La  orbitals  since  their  effect  is  only  evident  high 
above  E^.  In  this  study  we  examine  the  effect  of  each  of 
the  above  orbitals  on  how  accurately  the  first-principles 
band  structure  can  be  reproduced. 

In  this  work  we  fit  the  NRL-TB  method  to  an  LAPW 
band  structure  ranging  from  2.7  eV  below  to  0.8  eV  above 
ep-  The  TB  band  structure  fits  the  LAPW  results  very 
well  near  ep  and  perfectly  reproduces  the  Fermi  surface 
near  the  T  and  M  symmetry  points. 

The  TB  densities  of  states  are  also  in  very  good  agree¬ 
ment  with  the  corresponding  LAPW  results,  including 
a  comparison  by  orbital  decomposition.  In  addition,  we 
have  studied  the  variation  of  the  total  energy  with  re¬ 
spect  to  the  position  of  the  As  atoms.  We  have  found 
that  the  TB  total  energies  fit  the  LAPW  values  very  well 
even  for  energies  that  we  did  not  include  in  our  fit. 


II.  COMPUTATIONAL  DETAILS 

The  equilibrium  structure  of  non-magnetic  LaOFeAs 
is  the  AsCuSiZr  structure^.  This  structure  is  tetrag¬ 
onal,  space  group  P4/nmm,  with  eight  atoms  in  the 
unit  cell  and  lattice  constants  a  =  7.626  Bohr  and 
c  =  16.518  Bohr.  The  Oxygen  atoms  occupy  the  (2a) 
Wyckoff  positions,  and  Iron  the  (2b)  positions.  The  La 
and  As  atoms  occupy  (2c)  sites,  with  2  =  0.14154  for 
La  and  0.6512  for  As.  We  used  a  regular,  T-centered 
8x8x4  k-point  mesh,  which  results  in  225  points  in  the 
irreducible  part  of  the  Brillouin  zone.  The  LAPW  basis 


functions  were  cut  off  at  RKmax  =  8.5,  with  approx¬ 
imately  1250  basis  functions  at  each  k-point.  To  help 
convergence  we  broadened  the  spectrum  using  a  Fermi 
distribution  at  a  temperature  of  5  mRy. 

As  expected  our  LAPW  energy  bands  and  densities  of 
states  are  basically  identical  to  those  published  by  Singh 
and  Dr£.  We  also  performed  21  additional  LAPW  cal¬ 
culations  by  varying  the  above  structural  parameters  a 
and  for  the  purpose  of  creating  a  first-principles  total 
energy  database  to  use  in  our  TB  calculations. 

Our  TB  Hamiltonian  was  built  following  the  NRL-TB 
method.  We  summarize  below  the  basic  equations  of 
this  scheme  which  is  based  on  a  Slater-Koster  approach-^ 
with  two-center  parameters. 

Unlike  the  general  NRL-TB  method,  where  we  include 
s,  p ,  and  d  orbitals  for  each  atom,  here  our  basis  set 
includes  only  the  d  orbitals  for  Fe  and  the  p  orbitals  for 
As  and  O.  All  other  contributions  to  the  band  structure, 
including  the  effects  of  La,  are  ignored  as  they  have  little 
weight  in  the  region  between  the  bottom  of  the  As-0  p 
bands  and  until  well  above  the  Fermi  level.  Since  there 
are  two  Fe,  As,  and  O  in  each  unit  cell  of  the  structure, 
we  will  end  with  a  22x22  matrix  to  diagonalize  at  each 
k-point.  In  addition,  we  limit  ourselves  to  an  orthogonal 
Hamiltonian,  so  we  ignore  the  possible  overlap  hopping 
parameters. 

Onsite  Parameters  In  the  NRL-TB,  the  onsite  energies 
are  determined  by  the  interaction  of  an  individual  atom 
with  its  environment.  In  our  study  of  LaOFeAs,  however, 
we  are  considering  at  most  small  displacements  of  the 
atoms  around  the  equilibrium  positions.  We  therefore  use 
a  constant  value  for  the  onsite  parameters.  However,  we 
recognize  that  the  symmetry  of  the  iron  states  allows  the 
orbitals  to  have  different  onsite  parameters.  Accordingly, 
we  chose  four  onsite  parameters  for  Fe:  one  for  the  xz 
orbital,  which  will  be  equal  to  the  yz  parameter,  one  each 
for  the  xy,  x2  —  y2  and  3 z2  —  r2  orbitals.  We  also  have  one 
onsite  parameter  for  each  the  As  p  orbitals,  and  another 
for  the  O  p  orbitals,  for  a  total  of  six  independent  onsite 
parameters,  listed  in  Table  HI 

Hopping  Parameters  The  usual  spd  Slater-Koster 
scheme  has  10  two-center  hopping  parameters  for  like- 
atom  interactions  and  14  parameters  for  the  interaction 
of  two  unlike  atoms.  Since  we  limit  our  system  to  Fe  d 
and  As/O  p  orbitals,  we  have  only  the  following  interac¬ 
tions: 

Fe-Fe:  dda,  ddir,  ddd 

As- As  (or  O-O):  ppa,  ppn 

Fe-As  (or  Fe-O)  (Fe  atom  first):  dpa ,  dpir 

As-O:  ppcr,  ppir 

for  a  total  of  13  Slater-Koster  parameters.  As  with 
the  standard  NRL-TB,  these  parameters  depend  only 
on  the  distance,  R ,  between  the  two  atoms,  and  are 
parametrized  according  to  the  formula 

HU^(R)  =  (aep ii  +  bw^R  +  cu'^R2)  e^p(-d2e, ^R)F(R)  , 

(1) 
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TABLE  I:  Onsite  parameters  for  Fe-As-O,  determined  using 
the  methods  described  in  the  text.  All  energies  are  in  Ryd¬ 
bergs. 


Atom 

Orbital 

Onsite 

Fe 

yz  ( zx ) 

0.51108 

xy 

0.54617 

2  2 
x  -y 

0.54548 

3z2  -  r2 

0.5513 

As 

P 

0.18566 

O 

P 

0.39230 

where  a,  b ,  c  and  d  are  fitting  parameters,  and  the  cutoff 
function  F(R)  has  the  form 

F{R)  =  1/1  +  exp[(i?  —  R0)/£],R  <  Rq 

=  0,R>  R0  ■  (2) 

In  our  fits  we  take  R  =  14  Bohr  and  l  =  0.5. 

The  parameters  in  the  above  equations  are  determined 
by  a  least-squared  fit  to  the  LAPW  eigenvalues  at  225 
k-points  in  the  irreducible  Brillouin  zone.  We  have 
found  that  applying  group  theory  to  block-diagonalize 
the  Hamiltonian  at  high  symmetry  points  is  essential  for 
obtaining  a  good  fit. 

III.  RESULTS 

The  onsite  parameters  determined  from  our  fitting  pro¬ 
cedure  are  given  in  Table  [U  To  make  our  results  easier  to 
compare  to  first-principles  calculations  we  have  rotated 
the  Cartesian  axis  by  45°  relative  to  the  primitive  vectors 
of  the  tetragonal  unit  cell,  so  that  the  lobes  of  the  x2  —  y2 
Fe  orbitals  point  to  the  nearest-neighbor  Fe  atoms. 

In  Table  |n]  we  list  the  Slater-Koster  hopping  parame¬ 
ters  generated  by  Eq.  [T| 

These  parameters  are  designed  to  fit  the  9th  to  21st 
bands  of  our  TB  Hamiltonian  to  first  principles  results 
for  three  nearest  neighbors  of  Fe-Fe,  As- As,  O-O,  Fe-As, 
Fe-O,  and  As-0  interaction. In  TableHIDwe  show  the  RMS 
error  per  band  for  this  fit,  indicating  the  high  quality  of 
the  fit  to  the  LAPW  eigenvalues.  It  should  be  noted  that 
near  the  Fermi  level  which  is  in  the  vicinity  of  the  17th 
and  18th  bands  the  deviation  between  LAPW  and  TB 
is  on  the  average,  including  all  the  225  k-points,  3mRy. 
At  the  high  symmetry  points  such  as  at  P,  X  and  M  the 
deviation  is  less  than  ImRy,  as  can  be  seen  in  Fig.[l]which 
shows  a  comparison  between  LAPW  and  TB  results. 

The  aforementioned  13  bands  span  an  energy  range 
from  0.4Ry  to  0.75Ry  as  shown  in  Fig.  [1]  In  this 
figure  one  can  see  that  from  0.45  Ry  where  a  gap  is 
present  to  0.70  Ry  the  fit  is  good  and  around  the  Fermi 
level(0.605Ry)  the  agreement  between  TB  and  LAPW  is 
excellent;  the  energy  bands  clearly  show  holes  around 
the  center  of  the  Brillouin  zone  and  electron  pockets 
around  the  high  symmetry  point  M.  This  suggests  that 


TABLE  II:  Parameters  used  to  construct  the  Slater-Koster 
hopping  parameters  ©,  as  described  in  the  text.  Energy 
units  are  in  Rydbergs,  and  distances  are  in  Bohr. 

Fe  Fe  R  dda  ddn  ddS 

5.329  -0.02771  0.01001  0.00031 

7.626  0.00546  0.00029  0.00750 

_ 10.784  0.00364  -0.00500  0.00008 

As  As  R  ppa  ppn 

7.350  0.05880  0.08276 

7.626  0.06633  0.04262 

_ 10.784  0.01041  -0.05779 _ 

OO  R  ppa  ppn 

5.392  0.01885  -0.00783 

7.626  0.00939  -0.00534 

_ 10.784  0.00208  0.00085 _ 

Fe  As  R  dpa  dpn 

4.558  0.17916  0.00931 
8.884  -0.00751  -0.02974 

_ 11.708  -0.00073  -0.00090 _ 

Fe  O  R  dpa  dpn 

8.259  -0.00319  -0.00338 
9.863  0.00021  0.00240 

_ 11.241  -0.01449  0.00648 _ 

As  O  R  ppa  ppn 

6.909  0.00513  -0.02238 
10.290  0.01562  -0.00206 
11.412  0.00591  -0.00028 


FeAsO:  TB  /LAPW  Band  comparison  (13  bands  fit) 


FIG.  1:  A  comparison  of  the  thirteen  LAPW  and  TB  bands. 


the  TB  Fermi  surface  preserves  all  the  characteristic  fea¬ 
tures  found  in  the  first-principles  calculations.  In  figures 
[2][5]we  show  a  comparison  between  TB  and  LAPW  densi¬ 
ties  of  states.  In  Fig.[2]the  total  DOS  are  compared  where 
one  can  see  that  TB  matches  well  the  two  LAPW  peaks 
below  £p  and  the  one  peak  above  £f ■  The  site  decom¬ 
posed  DOS  shown  in  Fig.  [2]  shows  very  good  agreement 
for  the  Fe  d-states,  and  reasonable  to  good  agreement 
for  the  As  and  O  p-states.  One  can  conclude  that  the 
TB  produces  reliable  results  not  only  for  the  eigenvalue 


TABLE  III:  RMS  ERROR:  per  band  for  225  k-points 
Band  RMS  Error  (Ry) 
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9 

0.008365 

10 

0.009243 

11 

0.011206 

12 

0.009935 

13 

0.011851 

14 

0.007552 

15 

0.005671 

16 

0.005491 

17 

0.002917 

18 

0.002912 

19 

0.007435 

20 

0.008658 

21 

0.012827 

22 

0.010253 

FIG.  2:  (color  online)  Comparison  of  the  LAPW  and  TB 
DOS.  Site  decomposition  for  Fe,  As,  and  O  is  also  shown. 

spectrum  but  also  for  the  eigenvectors. 

The  orbitally  decomposed  DOS  are  also  in  at  least 
semi-quantitative  agreement  between  LAPW  and  TB  re¬ 
sults.  Note  that  since  the  LAPW  angular-momentum 
components  of  the  DOS  are  projections  inside  the  muffin- 
tin  spheres,  an  exact  comparison  cannot  be  made.  Still, 
the  TB  is  a  powerful  tool  in  explaining  features  of  the 
DOS  that  cannot  be  addressed  by  the  first-principles  cal¬ 
culations.  For  example,  we  can  trace  down  the  origin  of 
the  the  so-called  pseudogap  above  Ef  shown  in  Fig.  [2] 
Since  the  local  environment  of  the  Fe  atom  is  tetragonal 
(distorted  to  some  degree  in  most  1111  type  compounds), 
and  given  that  the  strongest  hopping  parameter  is  the  Fe- 
As  dpa ,  as  seen  in  Table  m  one  might  expect  to  see  a 
gap  or  pseudogap  between  a  lower  eg  doublet  and  an  up¬ 
per  t2g  triplet,  commensurate  with  the  local  symmetry. 


FIG.  3:  (color  online)  The  DOS  with  only  nearest  neighbor 
Fe-As  hoppings  included.  The  expected  ligand  field  splitting 
corresponding  to  the  tetragonal  environment  of  the  Fe  atom 
is  clearly  visible.  As  a  comparison,  the  DOS  including  all 
hoppings  is  shown. 


However,  in  reality  the  gap  occurs  between  a  lower  peak 
in  the  DOS  containing  three  states  per  Fe  (six  total)  and 
and  upper  peak  containing  two  states  per  Fe  (four  total), 
i.e.  the  reverse  of  the  naive  ligand  field  (or  crystal  field) 
expectation  (this  can  be  seen  in  Fig.  [2]  in  the  Fe  part  of 
the  spectrum).  Here  we  use  our  TB  model  to  eliminate 
hoppings  one  by  one  and  trace  down  the  origin  of  the 
pseudogap. 

First  we  turn  off  all  but  the  Fe-As  nearest  neighbor 
hoppings  and  indeed  find  a  lower  doublet  and  upper 
triplet,  as  must  be  the  case  (see  Fig.  H.  The  splitting 
between  triplet  states,  and  to  a  lesser  degree  the  doublet 
states,  is  due  to  the  imperfect  tetrahedron. 

Next,  we  eliminate  all  Fe-As  hopping,  but  allow  hop¬ 
ping  between  the  two  crystallographically  inequivalent 
Fe  sites.  As  can  be  seen  in  Fig.  [4]  this  has  the  effect 
of  creating  bonding/anti-bonding  complexes  within  each 
orbital  designation.  Interestingly,  the  xz  +  yz  orbitals 
do  not  split  at  all,  remaining  as  a  single  peak,  while  the 
most  strongly  split  orbital  is  xy.  Already,  something  like 
a  pseudogap  can  be  seen  to  be  forming  around  0.57  Ry, 
with  a  spectral  distribution  quite  different  from  the  lig¬ 
and  field  gap  distribution  in  Fig.  [3l 

Finally,  we  allow  all  Fe-As  hoppings,  but  remove  direct 
Fe-Fe  hopping.  This  scenario  accounts  for  interaction  be¬ 
tween  the  two  Fe  sites  in  the  unit  cell  via  As.  In  Fig.  [5j 
a  bonding/anti-bonding  splitting  again  occurs  in  most  of 
the  orbitals  with  the  strength  of  the  splitting  consider¬ 
ably  stronger  than  that  initiated  by  Fe-Fe  direct  hopping. 
In  this  case  the  xy  +  xz  orbital  does  undergo  splitting, 
but  the  z2  and  xy  orbitals  do  not.  A  fairly  strong  pseudo¬ 
gap  emerges  with  the  majority  of  the  weight  in  the  lower 
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FIG.  4:  (color  online)  The  DOS  with  all  hoppings  eliminated 
except  direct  hopping  between  Fe  atoms.  Most  of  the  orbitals 
split  into  bonding/anti-bonding  combinations  with  only  the 
xz  +  yz  combination  remaining  as  a  single  peak. 


FIG.  5:  (color  online)  The  DOS  including  all  Fe-As  hoppings, 
but  without  any  direct  Fe-Fe  hopping.  The  z2  and  xy  orbitals 
remain  unsplit,  while  all  other  orbitals  show  rather  strong 
bonding/anti-bonding  peak  separations. 


complex. 

Combining  the  information  from  all  three  reduced  hop¬ 
ping  diagrams,  we  can  understand  how  the  pseudogap 
forms.  It  is  a  result  of  a  combination  of  strong  and  weak 
bonding/anti-bonding  splitting  of  the  orbitals,  due  to 
crystallographically  inequivalent  Fe  atoms.  Hopping  oc¬ 
curs  both  via  As  and  directly  between  Fe  atoms  (further 
hoppings  also  surely  contribute  somewhat).  In  Fig.[Gl  we 
show  how  very  strong  splitting  of  the  f2g-triplet  derived 
states  coupled  with  weaker  splitting  of  the  eg  doublet- 
derived  states  results  in  the  calculated  pseudogap.  The 
xz  +  yz ,  x2  —  y2  orbitals  split  strongly  into  two  peaks  per 


FIG.  6:  (color  online)  A  schematic  illustrating  the  origin  of 
the  pseudogap.  Each  of  the  original  five  orbitals  is  doubly 
degenerate  (one  set  coming  from  each  Fe  in  the  unit  cell). 
The  xz  +  yz,  x2  —  y2,  and  xy  orbitals  bonding  and  anti-bond 
with  sufficient  strength  to  form  upper  and  lower  complexes. 
The  z2  orbital  has  a  small  splitting,  but  both  peaks  remain 
in  the  lower  complex.  The  pseudogap  is  between  the  upper 
(anti-bonding  only)  and  lower  (most  anti-bonding)  complexes. 


orbital,  forming  an  upper  complex  of  three  states  and  a 
lower  one  also  of  three  states.  The  xy  orbital  has  a  weaker 
splitting,  but  strong  enough  still  to  place  one  state  in  the 
upper  complex  and  one  in  the  lower.  Not  surprisingly,  the 
z2  orbital,  which  is  pointed  mainly  out-of-plane,  has  the 
weakest  splitting  such  that  both  the  bonding  and  anti¬ 
bonding  peaks  remain  in  the  lower  complex.  Thus,  there 
are  six  states  in  the  lower  complex  and  four  in  the  upper. 
Note  that  the  actual  splitting  of  the  xy  state  is  somewhat 
more  complicated  than  the  simple  bonding/anti-bonding 
schematic  suggests.  This  simply  mirrors  the  fact  that 
the  DOS  itself  is  not  composed  of  two  simple  peaks  of 
precisely  six  and  four  states  each,  but  rather  has  some 
secondary  peak  structures.  By  turning  off  specific  hop¬ 
pings,  we  have  elucidated  from  an  atomic  orbital  point 
of  view,  the  mechanism  that  gives  rise  to  the  pseudogap. 
It  requires  both  direct  Fe-Fe  interaction  and  interaction 
via  the  intermediate  As. 

The  second  fit  that  we  performed  has  the  objective  to 
fit  the  volume  and  the  As  position  variations  of  the  total 
energy.  The  height  of  the  As  ion  has  strong  effects  on  the 
electronic  structure  and  magnetism  and  may  even  be  able 
to  switch  the  pairing  symmetry7!17!18.  For  this  purpose 
we  run  21  separate  LAPW  calculations  to  fit  in  our  TB 
scheme  the  total  energy.  In  this  calculation  we  fitted  all 
the  22  energy  bands  that  correspond  to  our  Fe(d)-As(p)- 
O(p)  TB  Hamiltonian.  This  fit  does  not  give  the  excellent 
fit  to  the  LAPW  Fermi  surface  that  we  found  in  the  first 
fit.  However,  it  fits  the  21  LAPW  total  energies  perfectly 
with  an  RMS  error  of  0.0004  Ry.  In  Fig.  [7]  we  show  the 
total  energy  results  of  this  fit.  We  used  these  TB  pa¬ 
rameters  to  calculate  total  energies  outside  our  database 
and  compared  with  independent  LAPW  results  not  in¬ 
cluded  in  the  fit  and  found  very  good  agreement.  Using 
these  results  we  found  the  lattice  equilibrium  parameters 
to  be  a=7.4  a.u.,  c=16.4  a.u  and  As  position  z=0.64,  in 
reasonable  agreement  with  the  experimental  values. 


6 


a  (Bohr) 

FIG.  7:  (color  online)  The  total  energy  curves  generated  by 
TB  fits  (open  symbols)  to  LAPW  data  (solid  symbols)  for 
various  values  of  the  As  height  above  the  Fe  plane  and  planar 
lattice  constant,  for  fixed  lattice  constant  c  =  16.4  a.u. 


IV.  CONCLUSIONS 

We  report  TB  results  on  LaOFeAs  obtained  by  the 
NRL-TB  method  via  a  fit  to  LAPW  eigenvalues  and  to¬ 


tal  energies.  Two  TB  parametrizations  were  performed: 
the  first  aims  at  reproducing  the  energy  bands  in  an  en¬ 
ergy  range  from  2.7  eV  below  to  0.8  eV  above  Ef  with 
superior  accuracy  around  Ef  ■  From  this  parametrization 
an  analysis  of  the  orbital-decomposed  DOS  shows  that 
the  mechanism  which  creates  the  pseudogap  above  Ef 
comes  from  a  direct  Fe-Fe  interaction  and  from  an  Fe-Fe 
interaction  through  the  intermediate  As  atom.  In  our 
second  TB  parametrization  we  focus  on  the  energetics 
of  the  LaOFeAs  system  finding  a  TB  Hamiltonian  that 
fits  the  LAPW  total  energies  as  a  function  of  volume  and 
As  position  very  well.  We  propose  that  this  TB  Hamilto¬ 
nian  will  be  very  useful  in  carrying  out  many-body  theory 
with  a  more  realistic  Hamiltonian  than  those  employed 
previously  containing  just  the  d-iron  orbitals. 
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